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Abstract 

We report on a high statistics quenched lattice QCD calculation of the deep-inelastic 
structure functions Fi, F2, gi and ^2 of the proton and neutron. The theoretical basis 
for the calculation is the operator product expansion. We consider the moments of 
the leading twist operators up to spin four. Using Wilson fermions the calculation 
is done for three values of k, and we perform the extrapolation to the chiral limit. 
The renormalization constants, which lead us from lattice to continuum operators, are 
calculated in perturbation theory to one loop order. 
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1 Introduction 



Deep-inelastic lepton-nucleon scattering is described by four structure functions: Fi, F2, 
gi and §2- The spin- averaged structure functions Fi and F2 carry information about the 
overall density of quarks and gluons in the nucleon. They have played a seminal role in 
the development of our current understanding of the structure of hadrons. The polarized 
structure function gi goes one step further and probes the distribution of quarks of a given 
helicity in the longitudinally polarized nucleon. Recent measurements of gi [1] have revealed 
the (at first sight) surprising result that only a small fraction of the nucleon's spin is carried 
by quarks. This has triggered a great deal of interest in the subject. The other polarized 
structure function g2 has no interpretation in purely partonic language. It involves a twist- 
three operator and thus offers the first direct measurement of higher twist operator matrix 
elements [2]. Experiments that measure g2 are currently being performed at DESY and 
SLAC. 

We have initiated a program to compute Fi, F2, di and g2 on the lattice [3]. Eor an 
earlier attempt to compute the unpolarized structure functions see [4]. The theoretical basis 
for such a calculation is the operator product expansion (OPE), which relates the moments 
of the structure functions to forward nucleon matrix elements of certain local operators. 
Where a parton model interpretation exists, it can be mapped onto an OPE analysis. Our 
calculation will be in the quenched approximation, where internal quark loops are neglected. 
In this paper we shall also neglect gluonic operators, which contribute only to higher order 
in the coupling constant expansion. Eor the unpolarized structure functions we then have 
for the leading twist contribution 

2 j' dxx--'F,{x,Q') = J2 ^i,lif'VQ\9{f^))vl/\fi), 
j'dxx-'F2{x,Q') = J: i%'/Q'^9{p^))^i^\l^^) (1) 



f=u,d 



for n > 2 (generally even n), where 
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^J2iP^^\^{tli-^n}\P^^ = '^^n ^iPi^i ■■■Pi^n- traces]. 



^l^U = y V'7..^..---^.„V'- traces (2) 

with t/j = u[d) for / = m(J), and 

c[%'/Q\9{p^)) = Q^^^'{l+g{jifcrA/IQ\g{m. 
ci%VQ\9if^)) = Q^^^\l+9{}ifc2A/IQ\9{m- (3) 

Here // denotes the subtraction point, and {• • •} indicates symmetrization. We have chosen 
the normalization {p, s\p' , s') = {2'KY2Fj^6{p — p')6g^gi, = — m^. The moments of Fi, F2 
have the parton model interpretation 

^(/) = (4) 
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where x is the fraction of the nucleon momentum carried by the quarks. In the quenched 
approximation the above-mentioned equations hold for odd n as weU. 

For the polarized structure functions we have, again for the leading twist contribution, 

2 I" dxx-g,{x,Q') = Ij2 e[%'/Q\g{fi))al!\fi), (5) 
Jo z , , 

j=u,d 

2 C dxx-g,{x,Q') = \^ Y: [4%VQ\9if^)) " e[%VQ\g{fi)) a(^)(/.)] (6) 

Jo 2 n + 1 f , 

f=u,d 

for even n and n > (n > 2) for gi (5'2), where [2] 

?\ ^ ^ nU)\ J, 



{P,^\^{at!^-^r.}\P^^ = ~^^n^^^P^^ ■■■P^n + traces] 

{P^^\^[1{L]---^„}\P^^ = -^^-p^dl^^ii^^Pf^i - ■5mP^)Pa*2 ■■■Pf^n + traces], (7) 

O'Jil,^ = (0''V7.75^..---^.„V'- traces 
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e[%VQ\9{^)) = Q^^^'{l+g{jiferA/IQ\9{m. 
ei%VQ\9if^)) = Q^^^\l+9{}ife2A/IQ\9{m- (8) 
In eq. (7) [• • •] indicates antisymmetrization. In parton model language 

aj") = 2Au, a^a^ = 2Ad, (9) 

where Am, Ad determine the fraction of the nucleon spin that is carried by the quarks. 
A similar interpretation holds for the higher spin operators. The structure function g2 
consists of two contributions: a^^"^ is the so-called Wandzura-Wilczek contribution [5] which 
corresponds to twist two, whereas d^/'^ corresponds to twist three. In addition to d^/'^ there 
is a contribution to g2 from a twist-three operator which is proportional to the quark mass 
[6]. Since we are mainly interested in the chiral limit we have neglected that. 

2 Lattice Calculation 

We perform our quenched QCD calculations for Wilson fermions with r = 1 on a 16"^ X 32 
lattice. For the gauge coupling we take /3 = 6/g'^ = 6.0. To be able to extrapolate our 
results to the chiral limit, we run at three different hopping parameters, k, = 0.155,0.153 
and 0.1515. This corresponds to physical quark masses niq of roughly 70, 130 and 190 MeV, 
respectively, using the perturbative relation 

m,a = 0.56(- - — ), (10) 

where a is the lattice spacing. 
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For the gauge update we use a cycle consisting of a single three-hit Metropolis sweep 
followed by 16 overrelaxation sweeps [7]. We repeat this cycle 50 times in order to generate 
a new configuration. The resulting configurations seem to be independent. We see no 
correlations between hadronic quantities calculated on different gauge field configurations. 
The calculations are carried out on Quadrics Q16 and QH2 parallel computers. For details 
of the implementation of our code on these machines see [3]. So far we have collected of the 
order 1000, 600 and 400 independent configurations at the three n values. 

To calculate the nucleon matrix elements we first compute two- and three-point correla- 
tion functions defined by 

a,fj 

Cr{t,T,p,0) = J2^rj,4B4t,p)O{T)Brj{0,p)). (11) 

a,f3 

The lattice operators O used here are obtained from the operators in the euclidean contin- 
uum, which look exactly like the operators in eqs. (2), (7) up to factors of i, by replacing 
the covariant derivative by the lattice covariant derivative so that D ^ becomes 

Dt,{x,y) = ^[U^,{x)Sy^^+n - ujix - fi)Sy,^-n]. (12) 

We write the ratio of three- to two-point correlation functions as 

i?(t,T,p,r,(9) = Cr(t,T,p,(9)/Ci(^^^^)(t,p) 

= ^Tr^^F(r,J) (13) 

(for < T < t) with 

F{V,J) = U:v[VNJN], 

N = j4-iP7/Ep-+mN/Ep- (14) 

and J7 defined by 

(p, s\0\p, s) = u{p, s)Ju{p, s). (15) 

When calculating three-point functions it is particularly important that the baryon op- 
erator B has only little overlap with excited baryon states, in order to make the plateau 
region in r as broad as possible. As our basic proton operator we use (with C = 7472 in our 
representation) 

Bait.P) = E e-'P"e,,,ul{x){u'{x)Cjsd'{xj) (16) 

x,a,b,c 

with two important improvements. First we use 'Jacobi smearing' [8] (a version of [9]) in 
order to have an extended proton operator. Thus each quark operator in eq. (11) is replaced 
by 

JV, _^ 

i^^i^' = Y.i^sDri^, (17) 

n=0 
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Figure 1: (a) Effective nucieon mass piot for p = (bottom), and effective nucieon energy 
piot for \p\ = 27r/f6 (top) at k, = 0.f55. Botli source and sink; are smeared. Tlie liorizontai 
fines indicate the resuit of the fit as weii as the fit intervai. (b) m^, and function 
of 1/k, together with other recent resuits of the iiterature: Q this work;, □ Ref. [fO], O 
Ref. [f f], A Ref. [12]. The fit to aii data points combined gives Kc = O.f 5699(5). 
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0.1515 


n 

0.153 


0.155 


nip 
rriN 


0.504(2) 
0.570(2) 
0.900(5) 


0.422(2) 
0.507(2) 
0.798(5) 


0.297(2) 
0.422(2) 
0.658(5) 



Table 1: The hadron masses in lattice units at /3 = 6.0. 

and similarly for as we smear both source and sink. We found suitable values of the 
parameters to be Ng = 50, = 0.21, which for our largest n value gave a rms radius of 
about 4, corresponding roughly to 0.5 fm, i.e. half the nucleon radius. Secondly we replace 
each spinor by 

^ V"^^ = \{i + 74) V-, V ^ V''^'' = 1(1 + 74). (18) 

This replacement leaves quantum numbers unchanged, but we expect it to improve overlap 
with those baryons which have slow-moving valence quarks. Practically this means that for 
each baryon propagator we invert on a smeared local source and consider only the first two 
Dirac components. So we only have 2x3 inversions to perform rather than the usual 4x3 
inversions, which saves a factor two in computer time in the inversion. 

In Ref. [3] we have seen that the projection (18) is particularly effective at reducing 
unwanted backward propagating states, which extends the window that one can practically 
use for matrix element calculations to well above half the temporal extent of the lattice. In 
Fig. la we show a plot of the resulting effective nucleon energy, as given by ln(C(t)/C(t-|-l)), 
for p = and \p\ = 27r/16 at our smallest quark mass. For zero nucleon momentum we find 
a good plateau with a proton mass of 0.658(5). For the lowest non-zero momentum we find 
an energy of 0.765(11), which is in good agreement with the continuum dispersion relation. 
In both cases we see that after a distance of about four time units there is very little trace 
of an excited state. In Table 1 we give the mass values of the nucleon together with those of 
the TT and p for our three values of n. The chiral limit is obtained by extrapolating in Ij n 
to zero TT mass. Assuming, as usual, that depends linearly on 1/k, we obtain from our 
data the critical value Hc = 0.15693(4). In Fig. lb we plot m^, and function of 

Ij n. In this plot we have also included other recent results at /3 = 6.0. A combined fit gives 
m]\j jnip = 1.37 in the chiral limit. 

To calculate three-point functions we require additional propagators, one for each chosen 
t, p and r. We have fixed t at 13 and have chosen F = i(l-|-74), corresponding to the 
unpolarized case, and F = |(1 -|- 74)27572, corresponding to polarization (-|- - — ) in the 2- 
direction. For the momentum we have taken p = and p = (27r/16,0,0) = (pi,0,0). We 
have also considered t/j = d separately. This means that we must find 2x2x2 = 8 (half) 
quark propagators. The choice t = 13 is sufficient. Larger values of t lead to unacceptably 
large errors in the signal for R. Test runs for t = 17 turned out to have 0(2) larger errors, 
which roughly corresponds to the increase in the noise in the baryon correlation function 
from t = 13 to t = 17. 
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Representation 


















Ket. [ioj 


Ref. [14] 


c 


V2,a 


C{14} 


(6) 
^3 


6(+) 


+ 




^3^ 










V2,b 


C{44} - UO{ll} + 0{22} + 0{33}) 




1 




+ 


V3 


C{114} — |(C{224} + C'{334}) 


(8) 


8(+) 


— 






r? 






(+) 




V4 


C{1144} + C'{2233} — C'{1133} " C'{2244} 










+ 


ao 


01 


M 


(i 


1 
2 




+ 




•^{214} 


^3 


(i 


1 
2 




+ 


d2 




(8) 


8(+) 


+ 



Table 2: The lattice operators and their representation. The momentum is taken to be 
p = (27r/16, 0, 0) = (pi, 0, 0) in the case of U2,a, ^^3, ^^4, ^2, (^2 and p = elsewhere. C denotes 
charge conjugation. 

3 Lattice Operators and their Renormalization 

The bare lattice operators, (9(a), are in general divergent. We define finite operators 
renormalized at the scale by 

0{^i) = Zo{{a^i)\g{a))0{a), (19) 

where 

{q{p)\0{iJ:Mp)) = {q{p)\0{aMp)) (20) 

with \q{p)) being a quark state of momentum p. In the limit a — )• this definition amounts 
to the continuum, momentum subtraction renormalization scheme. 

The lattice operators transform under the discrete hypercubic group -ff(4) [13, 14]. They 
must be constructed such that they belong to a definite irreducible representation of the 
latter. In particular they must not mix with lower-dimensional operators. This is prerequisite 
to the operators being multiplicatively renormalizable. Furthermore, from the more practical 
point of view, the operators should only require a non-zero spatial momentum in at most 
one direction. We have considered the operators listed in Table 2. For the group theoretical 
classification of the lattice operators see Ref. [15]. The calculation of U2,a, ^^3, ^^4, 0.2 and 
^2 requires non- vanishing nucleon momenta. Note that for the quenched theory there is 
no mixing with gluon operators. In the continuum limit the matrix elements V2^a and ^2,5 
should be equal. At finite lattice spacing this provides us with a consistency check and gives 
information about possible lattice artifacts. 

We have computed the renormalization constants for our operators in the quenched 
approximation for Wilson fermions in perturbation theory to one loop order. For this task 
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(O) 


10 


Bq 


^0(1,5 = 1-0) 


rM5 
^0 


V2,a 
V2,b 

V3 {}{}: 
{}(): 

V4 
ao 
a.2 
d2 


16 

3 

16 

3 
25 
3 

U 

157 
15 



25 
3 
7 
3 


-3.165(6) 
-1.892(6) 
-19.572(10) 

0.3 ( 0(10j 

-37.16(30) 
15.795(3) 
-19.560(10) 
-15.680(10) 


1.0267(1) 
1.0160(1) 
1.1653(1) 

-0.0031(1 j 

1.314(3) 
0.8666(0) 
1.1652(1) 
1.1324(1) 


40 

9 

40 

9 

67 
9 

2216 
225 



67 
9 

13 
12 



Table 3: The renormalization constants in the quenched approximation. The errors quoted 
are a conservative estimate of the uncertainties in the numerical evaluation of the integrals 
involved. The numbers in the rightmost column represent the contribution of the continuum 
operators computed in the MS scheme. 

we have developed packages of computer algebraic programs using Mathematica and Maple 
to such a level that all what is needed as input is to state the Feynman rules in symbolic 
form, both for the continuum and the lattice part of the calculation. We will summarize our 
results here. A detailed account of our calculation will be given elsewhere [16]. 

In the case of it turns out that the operator (9{ii4}. — |(C{224} + ^'{334}) mixes with 
the operator [15] 

^(114) ~ ■^{^{224) + C'(334)), = C^^j. + 0^1,^ — 20i,nn (21) 

under renormalization. This operator is of mixed symmetry, is traceless and corresponds to 
the representation 8^"*"^ C = — as well. Thus we have 

^{}(/^) = %{}^{}(«) + %()^()(«). (22) 

where we have used a short-hand notation for the operators in Table 2 and eq. (21). We 
write {Cf = 4/3) 

Zo{{a.iJ.Y,g) = I - -f^CFboHaiJ.) + Bo]. (23) 

This is to be interpreted as a matrix equation in the case of V3. Our results for the anomalous 
dimensions 70 and the Bq^s are given in Table 3 for r = 1. The renormalization constants 
^ and have been given before in the literature [17, 18, 19]. We agree with the results 
of these authors. In the case of V3 the off-diagonal component of Z is negligibly small. 
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The structure functions do not depend on but (x" ^) and Am, Ad do. In the following 
we shall quote our results for 

= Q^ = ^ 2GeV', (24) 

which eliminates the logarithms in the Wilson coefficients and renormalization constants, 
and we will denote Zoil^g = f -0) by Zq- The corresponding numerical values are also listed 
in Table 3. As the Wilson coefficients are generally computed in the MS regularization 
scheme, one needs to know the renormalization constants in this scheme too. In Table 3 
we state the contribution of the continuum operators computed in the MS scheme. The 
difference of the Bq^s then gives the result in the MS scheme. 

The renormalization constants receive contributions from five different types of diagrams: 
the vertex, the leg self-energy, the leg tadpole, the operator tadpole and the operator comb 
diagrams. The tadpole diagrams give by far the largest contribution to the renormalization 
constants. The leg tadpole contribution is the same for all operators. The operator tadpole 
contribution is proportional to the number of covariant derivatives and has opposite sign to 
the leg tadpole. Leg and operator tadpole diagrams cancel each other in V2. This accounts 
for the small values of By^ ^, By^ ^. In ao only the leg tadpole contributes. In all other cases 
it is the operator tadpole diagram which dominates. 



4 Structure Function Results 

The next step is to calculate the ratio of three- to two-point correlation functions i?, as given 
in eq. (13), for the operators listed in Table 2. To make sure that we are computing the 
matrix elements of the lowest-lying state, i.e. the nucleon, we must look for plateaus in r, 
the time distance of the operator from the source, for ^ r ^ t = 13. In Fig. 2 we show R 
as a function of r for six of our operators at n = 0.153. The ratio Ry^ ^ not shown here is of 
the same quality as Ry^ ^. We find in all cases that the signal is practically constant for time 
distances larger than two lattice spacings from the source and from the sink. For 13 < r the 
signal is practically zero as one would expect. The fit interval is taken to be 4 < r < 9. The 
result of the fit is shown by the horizontal lines, and the errors are indicated by the dotted 
lines. 

The renormalized operator matrix elements are obtained from the ratio R by 

if II 

Rv2,a = 7^Pl^^2,a, Rv2,t = 7^mjvU2,6, 



1 1 o „ 11 

(25) 



1 1 mjv „ 111 

2,K ZBp^ Za^ 2k 6 
111 

Zd^ ZK 3 
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Figure 2: The ratio R for u and d quark insertions for k, = 0.153. (a) -'iRv2a^ (t*) Rv^^ (c) 
Rv^, (d) -ii?ao7 (e) -Ra2 s-iid (f) i?d2- Both source and sink are smeared. The source is at 
t = 0, the sink at t = 13. 
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We have defined the continuum quark fields by \/2/€ times the lattice quark fields. For the 
renormalization constants we take the perturbative values given in Table 3 and Z^^ = 1. 
In the case of V3 we have also computed the nucleon matrix element of the operator in 
eq. (21). It turned out to be noisy and consistent with zero within an error of roughly 
1/5 the magnitude of the leading symmetric contribution. Given the small off-diagonal 
component of the renormalization constant, we may thus safely neglect the effect of mixing. 
Tadpole resummation [20] would leave Zy^^^ ^V2b unchanged, while it would change the 
other renormalization constants by a few percent. The exact amount depends on how it 
is implemented, and there is considerable freedom to do so. (Better is to compute the 
renormalization constants non-perturbatively [21] what we are doing now [22].) The results 
are plotted in Fig. 3, and the numerical values are listed in Table 4. All our results are given 
for the proton. The distribution functions of the neutron are obtained by interchanging u 
and d. 

We shall now discuss our results in detail. The first important observation to make is 
that the values of (a;)^ and (x);,, which are obtained from different representations of the 
hypercubic group -ff(4) (cf. Table 2), are consistent with each other, within the error bars. 
This indicates that lattice artifacts are presumably small. A second observation is that all 
matrix elements show roughly a linear behavior in 1/k, i.e. in the quark mass (cf. eq. (10)). 
The lines shown in Fig. 3 are linear fits to the data. The result of the extrapolation is 
indicated by the solid circles and the solid box, and the numerical values of the fit are given 
in the last column of Table 4. 

Let us concentrate on the moments of the unpolarized structure functions first. We 
see that the lowest moment (n = 2) is practically independent of the quark mass. For 
growing n the moments show a stronger and stronger increase with the quark mass. For the 
distribution function itself this means that at small x quark mass effects are negligible, while 
at intermediate and large x its shape depends strongly on the magnitude of the quark mass. 
In the limit of large quark masses the higher moments approach the predictions of the non- 
relativistic quark model. In particular we find 2(a;"~^)('') for all n. In the chiral 

limit the picture changes completely. Whereas at small x the ratio of u to J distribution is 
roughly two, the ratio increases rapidly for larger values of x. 

We may compare our results in the chiral limit with the phenomenological valence quark 
distribution functions. In Fig. 3a-c we show the results of a recent such fit [23]. For n = 2 
the lattice values are significantly larger than the phenomenological values, for n = 3 they 
are about equal within the statistical errors, and for n = 4 the lattice values are significantly 
smaller than the phenomenological values. This holds for both, u and d quark distributions. 
Thus, our calculation predicts a valence quark distribution function that is more singular at 
small X than the phenomenological one, i.e. has a somewhat larger Regge intercept. At the 
moment we have no explanation for this discrepancy. According to our calculation 64% of 
the nucleon's momentum is carried by the (valence) quarks. 

We shall now turn to the discussion of the polarized structure functions. Let us first 
focus on Am and AJ, which in the quenched approximation determine the fraction of the 
proton spin that is carried by the valence quarks. Sea quark effects may be neglected for 
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6.35 6.40 6.45 6.50 6.55 6.60 6.65 

(a) 



6.35 6.40 6.45 6.50 6.55 6.60 6.65 

(b) 



-1 - 

6.35 6.40 6.45 6.50 6.55 6.60 6.65 
(C) 




Figure 3: The moments of the proton structure functions as a function of together with 
a hnear fit to the data. The sohd symbols indicate the extrapolation to the chiral limit. In (a) 
circles refer to (a;)^, boxes to (a;)^. In (a-c) we compare our results with the phenomenological 
valence quark distribution of Ref. [23] (fit -D-). The phenomenological moments are marked 
by asterisks. In (d) we compare our numbers with the phenomenological values of Ref. [24], 
which are marked by asterisks as well. (The value for Am is hidden behind the lattice 
number.) 
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Observable 


0.1515 


0.153 


K 

0.155 


/c - 1569 




0.477(23) 


0.436(20) 


0.457(28) 


0.430(43) 




0.478(22) 


0.473(20) 


0.475(20) 


0.473(32) 


\-^/av. 


0.478(16) 


0.455(14) 


0.466(16) 


0.452(26) 




0.226(11) 


0.201(10) 


0.195(12) 


0.174(20) 




0.226(10) 


0.219(09) 


0.211(09) 


0.204(15) 


\-^/av. 


0.226(08) 


0.210(07) 


0.203(07) 


0.189(12) 


\ / 


1 47n 1 
0.066(05) 


1 27n 0^1 
0.056(05) 


1 99(] 3^1 
0.047(06) 


1 04r20'l 
0.037(10) 




0.060(05) 
0.026(03) 


0.049(05) 
0.018(03) 


0.035(08) 
0.008(04) 


0.022(11) 
-0.001(07) 


Am 


0.938(45) 
-0.250(12) 


0.935(44) 
-0.250(12) 


0.863(43) 
-0.246(14) 


0.830(70) 
-0.244(22) 




0.170(10) 
-0.037(03) 


0.165(11) 
-0.037(04) 


0.161(20) 
-0.051(10) 


0.154(27) 
-0.050(12) 


4^) 


-0.110(05) 
0.020(02) 


-0.138(07) 
0.024(02) 


-0.201(16) 
0.036(06) 


-0.233(20) 
0.040(07) 



Table 4: Structure function results for the proton. All numbers refer to the momentum 
subtraction scheme. 

heavy quarks, and they drop out in the difference Am — AJ. In the chiral limit we obtain 

Am- AJ = 5^ = 1.07(9). (26) 

This is to be compared with the experimental value of the axial vector coupling constant 
qa = 1.26. In Fig. 3d we compare Am and A J individually with a recent phenomenological 
fit of the polarized structure function data [24], which naturally includes sea quark effects. 
If we add the sea quark contribution to our results - a recent lattice calculation [25] finds 
Am = A J = —0.14(5), As = —0.13(4) using perturbative renormalization factors - we 
would favor a somewhat smaller value for Am than the fitted value. For the total quark spin 
contribution to the nucleon spin we would furthermore obtain AE = 0.18(8), in agreement 
with the result of a full QCD calculation [26], i.e. including dynamical quarks. For heavy 
quark masses we find Am 1 and A J —1/4, in good agreement with the three-quark 
model [27]. 

By comparing the moments ao = 2A.q and a2 with those of the unpolarized structure 
functions we find that in the chiral limit gi is less singular than Fi as x goes to zero. This 
is also what one finds experimentally [28]. In the limit of large quark masses, on the other 
hand, it seems that gi is proportional to F-^. 
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If we combine our results with the perturbatively known [29] Wilson coefficients we can 
compute the moments of gi. In the chiral limit we obtain for the lowest moment 



Remember that 2GeV^. For the difference of proton and neutron structure functions 

we find 

/' dx{gl{x, Q') - g^{x, Q')) = 0.174(15). (28) 
Jo 

This is the quantity which should be compared with experiment, because here the sea quark 
contribution drops out. Our result is in good agreement with the phenomenological analysis 
[24, 28]. In the higher moments of gi sea quark effects should not play any role anymore 
either. In the chiral limit we obtain 



1 - - ' 0.0150(32) proton 

I -0.0012(20) neutron 



d..%M-)= ■ ; (29) 



Here we have converted the renormalization constants to the MS scheme, because the Wilson 
coefficients were computed in this scheme too. This result is consistent with experiment [30]. 

Let us finally discuss the structure function g2. From Fig. 3f we read off that the twist- 
three contribution is strongly mass dependent. While approaches zero in the heavy 
quark limit, for both u and d quark insertion, it is of the same order of magnitude as its 
twist-two counterpart 02 for small quark masses. In the chiral limit we obtain 



S,,2 J -0.0161(16) - 0.0100(22) = -0.0261(38) proton, 

^uxx(/2^x,v; I -0.0013(09) + 0.0009(13) = -0.0004(22) neutron. 



As before, Wilson coefficents [29] and renormalization constants are consistently computed 
in the MS scheme. In eq. (30) the first number comes from J2, while the second number 
comes from 02 (cf. eq. (6)). We see that the twist-three operator provides the dominant 
contribution. The Wandzura-Wilczek description of g2 [5] is a valid approximation for large 
quark masses, but for light quark masses it is definitely not. Our results seem to be in 
disagreement with recent estimates based on sum rules [31], which suggest that for the 
proton ^2 is very small. 



5 Conclusion 

We have presented results of a calculation of the lower moments of the polarized and unpo- 
larized deep-inelastic structure functions of the nucleon. The calculation has been performed 
in the quenched approximation, where sea quark effects are neglected, and it was done for 
three different quark masses. This allowed us to extrapolate our results to the chiral limit. 

The valence quark distributions that we have obtained differ somewhat from the phe- 
nomenological ones [23]. One explanation could be that at smaller values of higher twist 
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contributions are non-negligible, which have not been included in the phenomenological 
analysis. We plan to investigate this possibility in the future. Our results for the polarized 
structure functions are consistent with experiment, as far as data exist. A surprise was that 
the twist-three operator contributed so much to g2. 

It was interesting to see how the results varied with the quark mass. At large quark 
masses our results agree largely with what one would expect on the basis of the quark 
model. For small quark masses there are, however, significant changes. 

With the (raw) lattice data being relatively accurate now, the calculation of the renormal- 
ization constants has become a major issue. So far we have computed the renormalization 
constants in perturbation theory to one loop order. We hope to do better in the near fu- 
ture [22]. 

The renormalization constant for has independently been computed by the Rome group 
[32]. We have compared our results with theirs at intermediate stages of the calculation, and 
they agreed. These authors use a slightly different basis of operators from ours though. 
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